Charge disproportionation and collinear magnetic order 
in the frustrated triangular antiferromagnet AgNiC>2 
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We report a high-resolution neutron diffraction study of the crystal and magnetic structure of 
the orbitally-degenerate frustrated metallic magnet AgNiC>2. At high temperatures the structure is 
hexagonal with a single crystallographic Ni site, low-spin Ni 3+ with spin-1/2 and two- fold orbital de- 
generacy, arranged in an antiferromagnetic triangular lattice with frustrated spin and orbital order. 
A structural transition occurs upon cooling below 365 K to a tripled hexagonal unit cell containing 
three crystallographically-distinct Ni sites with expanded and contracted NiOe octahedra, naturally 
explained by spontaneous charge order on the Ni triangular layers. No Jahn- Teller distortions occur, 
suggesting that charge order occurs in order to lift the orbital degeneracy. Symmetry analysis of the 
inferred Ni charge order pattern and the observed oxygen displacement pattern suggests that the 
transition could be mediated by charge fluctuations at the Ni sites coupled to a soft oxygen optical 
phonon breathing mode. At low temperatures the electron-rich Ni sublattice (assigned to a valence 
close to Ni 2+ with S = 1) orders magnetically into a collinear stripe structure of ferromagnetic rows 
ordered antiferromagnetically in the triangular planes. We discuss the stability of this uncommon 
spin order pattern in the context of an easy-axis triangular antiferromagnet with additional weak 
second neighbor interactions and interlayer couplings. 

PACS numbers: 75.25.+Z, 71.45.Lr, 75.10.Jm, 75.40.Cx 



I. INTRODUCTION 

Two-dimensional frustrated quantum magnets have at- 
tracted a wide interest theoretically^ and experimentally 2 
as possible candidates to display strong fluctuations 
that could potentially stabilize unconventional ordered 
phases,— spin-liquid-, or orbital-liquid^ states. Low-spin 
triangular-lattice antifcrromagnets are canonical frus- 
trated models and the delafossite AgNiC>2 with Ni ions 
arranged in well-separated triangular lattices has been 
proposed to show frustration effects both in the mag- 
netic as well as in the orbital sector. Based on magnetic 
susceptibility measurements Shin et al£ proposed that 
Ni ions are in the low-spin configuration Ni 3+ (t2 g e g ) 
with 5=1/2, coupled by dominant in-plane frustrated 
antiferromagnetic interactions, and from x-ray measure- 
ments they proposed a high-symmetry crystal structure 
where each Ni ion has an unpaired electron in a doubly- 
degenerate e g orbital. The cooperative orbital order is 
strongly frustrated because the orbital interactions have 
a strong bond-directional dependence 7 favoring differ- 
ent orbitals for pairs of Ni ions along the three differ- 
ent directions in the triangular lattice and this leads to 
a large manifold of degenerate mean-field states. The 
ground state in this spin-orbital problem is still highly 
debated theoretically with proposals ranging from or- 
bital liquids to non-trivial forms of orbital order de- 
pending on fine details of the interactions^^ Among 



the experimentally explored candidate spin-1/2 materi- 
als to display this physics are NaTiC>2 (non-magnetic af- 
ter structural transition at low temperatures)^ NaNiC>2 
(quasi-2D spin-1/2 ferromagnet with ferro-distortive or- 
bital order and LiNiC>2 (no long-range magnetic or or- 
bital order, but evidence for local Jahn- Teller distortions, 
difficulty in preparing pure stoichiometric samples).- The 
triangular magnet AgNiO? and the two-silver-layer ver- 
sion Ag 2 Ni02 (Refs. [Illfl2T ) are relatively unexplored 
experimentally and promise to be rather different from 
the above-mentioned systems as they both show dom- 
inant antiferromagnetic interactions?^ thus potentially 
displaying frustrated magnetism. 

AgNiC>2 is part of the large family of delafossite 
A + M 3+ O2 materials often studied as candidate two- 
dimensional frustrated magnets because the transition 
metal ion (M) sits at the vertices of a triangular lat- 
tice in the basal plane, made up of a network of edge- 
sharing MOq octahedra. Like most delafossites, it oc- 
curs in two structural polytypes which differ in the way 
the MO2 layers are stacked along the c-axis: a 3-stage 
structure, where successive layers are in the same ori- 
entation but have an in-plane offset with Ni ions form- 
ing a 3-stage staircase along c-axis (3R polytype, rhom- 
bohcdral space group R3m, Refs. @Qj), or a 2-stage 
structure, where successive layers are stacked right on 
top of each other but are rotated by 180° (2H polytype, 
hexagonal space group P63 /mmc shown in Fig.QJa), Ref. 



Il5l ). All measurements reported here have been made 
on the less-studied 2-stage polytype, the so called 2H- 
AgNiC>2, recently synthesized using high-oxygen pressure 
techniques^ This polytype has metallic-like conductiv- 
ity from 300 K down to low temperatures and the sus- 
ceptibility indicates antiferromagnetic order near 20 Kji£ 
but the magnetic structure has not been determined up 
to now. 

In the ideal crystal structure of both 3R and 2H poly- 
types of AgNiC>2 Ni ions have both spin and orbital de- 
grees of freedom. The local crystal field is octahedral 
near-cubic and in the case of strong crystal field proposed 
herein&i£ the electronic state of Ni 3+ (3d 7 ) is the low- 
spin state t^eg with one unpaired electron (spin-1/2) in 
the upper e g level. A small trigonal distortion present in 
the crystal structure due to squashing of the NiC>6 octa- 
hedra along the c-axis changes the detailed wavefunctions 
of the orbital states but does not lift the two-fold degen- 
eracy of the upper e g level because it preserves a local 
three-fold symmetry rotation axis along c. Each Ni ion 
has a tendency to locally distort the environment to lower 
its orbital energy due to the Jahn- Teller effect, however 
the cooperative orbital order on the triangular lattice is 
frustrated as the orbital exchange favors occupation of 
different orbitals for pairs of Ni ions along the three dif- 
ferent in-plane directions. Such systems are susceptible 
to form an orbital liquid state at low temperatures or to 
have the orbital degeneracy lifted by structural distor- 
tions. 

In measurements reported here we find evidence for 
a weak structural modulation in 2H-AgNiC>2 leading to 
a tripling of the unit cell in the hexagonal basal plane. 
This can be naturally explained by charge disproportion- 
ation on the Ni sites into three sublattices, which we pro- 
pose occurs in order to lift the orbital degeneracy of the 
Ni 3+ ions. This physics is in sharp contrast to the insu- 
lator NaNiC>2 where the orbital degeneracy is lifted by 
Jahn- Teller orbital order leading to a monoclinic crys- 
tal structured We attribute this difference to the fact 
that 2H-AgNi02 being metallic charge transfer can be 
an energetically more favourable mechanism to lift the 
orbital degeneracy compared to local Jahn- Teller distor- 
tions found in more localized systems. 

At low temperatures the electron-rich Ni sublattice (at- 
tributed to a valence close to Ni 2+ with spin S = 1) 
orders magnetically in a collinear stripe structure with 
spins pointing along the c-axis and arranged in alter- 
nating ferromagnetic rows in the triangular plane. This 
magnetic structure cannot be explained at the mean-field 
level by a minimal spin model on a triangular lattice 
containing only nearest-neighbor antiferromagnetic ex- 
change and easy-axis anisotropy, and we propose that it 
is stabilized by additional weak second neighbor antifer- 
romagnetic in-plane interactions and/or weak ferromag- 
netic inter layer couplings. 

The rest of the paper is organized as follows. Next Sec- 
tion (JTTJ) presents the experimental details of the neutron, 
X-ray, susceptibility and specific heat measurements. 
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FIG. 1: (color online) Nominal crystal structure of 2H- 
AgNi02 deduced from X-ray measurements in Ref. [15 (space 
group P63/mmc-Dg h ). (a) There are two MO2 layers per unit 
cell related by a mirror plane reflection through the Ag + layer 
at z = 1/4. (b) Basal plane showing the triangular network of 
Ni ions (large red balls) coordinated by oxygens (small blue 
balls) . Thick solid line contour shows the unit cell and dashed 
line shows the unit cell tripling in the distorted structure. 



Diffraction measurements of the room-temperature crys- 
tal structure are presented and analyzed in Sec. IIIII where 
a lower-symmetry space group is proposed to accommo- 
date the observed structural modulation. A transition to 
the high-symmetry, undistorted crystal structure is ob- 
served upon heating to high temperatures and this is dis- 
cussed in Sec. IIVI followed by a symmetry analysis of the 
structural distortion in terms of symmetry-allowed basis 
vectors in Sec. |Vj This is used to propose that the mech- 
anism for the structural distortion is charge fluctuations 
at the Ni site coupled with a soft zone-boundary optical 
phonon involving oxygen breathing modes (further cal- 
culations using co-representation symmetry analysis to 
uniquely determine the distorted space group are pre- 
sented in Appendix B). Susceptibility and specific heat 
measurements are shown in Sec. IVI1 followed by mea- 
surements of the magnetic structure analyzed in terms of 
symmetry-allowed basis vectors and discussed in terms 
of a minimal Hamiltonian containing exchange and easy- 
axis anisotropy (Sec. IVII[) . Finally, the results are sum- 
marized and discussed in Sec. lVIlil For completeness we 
include in Appendix [5] a list of the measured supercell 
and magnetic structure factors. A partial account of the 
results describing the room-temperature crystal structure 
and low-temperature magnetic order has been reported 
in ref. fig. 



II. EXPERIMENTAL DETAILS 

Powder samples of the 2H-AgNi02 (< 1% admix- 
ture of the 3R polytype) were prepared from Ag20 and 
Ni(OH)2 using high oxygen pressures (130 MPa) as de- 
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TABLE I: Lattice parameters and atomic positions in the 
unit cell in the ideal (PG^/mmc) and the distorted (P6322) 
structural space groups at 300 K. Oxygen positions are pa- 
rameterized by an out-of-plane height zq — 0.08050(5) and an 
in-plane displacement e = 0.0133(2). The distorted unit cell 
is tripled in the hexagonal plane with an unchanged extent 
along the c-axis, but the origin is shifted by c/4 such that the 
two NiC>2 layers appear now at z = 1/4 and 3/4. Through- 
out this paper we use different symbols, ao and a — v3oo, 
to denote the hexagonal lattice parameter of the ideal and 
distorted structures, respectively. 
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scribed in Ref. UM Neutron diffraction patterns to 
probe the crystal and magnetic structure were collected 
using the high-resolution back-scattering time-of-flight 
diffractometers OSIRIS (0.65< Q <6 A" 1 ) and HRPD 
(2< Q <9 A^ 1 ) at the ISIS Facility of the Rutherford 
Appleton Laboratory in the UK. Preliminary measure- 
ments were also performed using GEM at ISIS and the 
monochromatic neutron diffractometer DIB at the In- 
stitute Laue-Langevin in France. The magnetic order 
parameter was obtained from elastic neutron scattering 
measurements using the direct-geometry time-of-flight 
spectrometer IN6 and the temperature-dependence of the 
lattice constants was also measured on D2B, both at the 
ILL. X-ray powder diffraction measurements to help solve 
the crystal structure at room temperature were made us- 
ing a Philips X'pert diffractometer (AcuK a = 1-54 A). 
Structural and magnetic refinement was made using the 
FullProf programme. 17 Susceptibility measurements were 
made using a SQUID magnetometer (Quantum Design 
MPMS) and specific heat data was collected on a pressed 
powder pellet using a Quantum Design PPMS system. 



III. CRYSTAL STRUCTURE 

The neutron powder diffraction pattern collected at 
room temperature (T=300 K, Fig. [2]) is overall in good 
agreement with the hexagonal space group PQ^/mmc 
proposed before. 15 However, a close inspection shows the 
presence of a number of additional low-intensity peaks 
(below 1% of the main peak) that could be indexed in 
this space group by fractional wavevectors such as (2/3,- 
1/3,0) and (2/3,-1/3,1) in Fig. O^b), and such super- 
cell peaks systematically accompany the main structural 



peaks throughout the wide Q-range probed [see panel 
a)] and are displaced in Q following the lattice contrac- 
tion upon cooling [see paired up-down panels b-c), e-f) 
and h-i)] . These extra peaks are naturally interpreted in 
terms of a structural modulation equivalent to a tripling 
of the unit cell in the basal plane. Complementary X- 
ray measurements [Fig. [3] d,g,j] did not show a mea- 
surable intensity at the supercell positions, suggesting 
that the structural modulation involves mainly displace- 
ments of the light oxygen ions which have a very small 
X-ray cross-section, as any significant displacements of 
the heavier Ag or Ni ions would have implied occurrence 
of supercell peaks not observed in the collected X-ray 
pattern. Therefore to model the distortion we refined 
only the oxygen positions and considered in order of de- 
creasing symmetry all subgroups of the ideal structure 
(P63 / mmc) compatible with a tripling of the unit cell in 
the ab plane, i.e. a unit cell of size V3ao x v3«o x c (6 
Ag, 6 Ni and 12 O atoms per unit cell). We eliminated 
the space groups that were not compatible with the ob- 
served low-temperature magnetic structure [see Sec. lVIIj 
where 2 Ni ions in the unit cell are magnetically ordered 
and 4 are unordered, i.e. we only considered the space 
groups where the magnetically-ordered and unordered Ni 
ions occupied distinct crystallographic sites. The high- 
est symmetry subgroup in which both the structural and 
magnetic data could be described is P6322 (no. 182) 
where all 12 oxygen atoms are symmetry-related and 
compared to the ideal structure arc displaced by a small 
amount e along one of the in-plane triangular directions. 
The observed supercell reflections could be well described 
by this model and the best fit to the data is shown in 
Figs. Hand 1 (solid lines, fl Bragg = 5.55%, R F = 6.49%) 
(structure factors are listed in Table. IVl in the Appendix 
A). The obtained lattice parameters and positions in the 
unit cell are listed in Table. I As a further test we re- 
fined the lattice parameters using only the main peaks 
or only the supercell peaks and obtained similar values 
(a = 5.09082(1) A, c = 12.24984(4) A and a = 5.0908(2) 
A, c = 12.250(2) A, respectively) corroborating the fact 
that the weak supercell peaks belong to the same phase 
as the main peaks and are not due to an extra phase. 

Fig. [4] shows a schematic diagram of the z = 1/4 
MO2 layer of the distorted structure. The distortion 
preserves the 3-fold rotation axis at each Ni site but the 
displacements of the oxygens lead to one Ni site with 
expanded Ni-0 bonds (Nil) and two other sites with 
contracted bonds (Ni2 and Ni3); the black circles cen- 
tered on the Ni sites correspond to Ni-0 distances of 
2.022 A (Nil) and 1.934 (Ni2 and Ni3), respectively. In 
Ni 2+ oxides 18 typical Ni-0 bond distances are about ~ 
2.09 A, whereas in Ni 4+ oxides^ they decrease to 1.92 
A. This comparison suggests a charge disproportionation 
among the Ni sites on the triangular layers in 2H-AgNi02 
between electron-rich Nil sites (expanded NiOg octahe- 
dron) with valence close to Ni 2+ and electron-depleted 
sites Ni2 and Ni3 (contracted octahedra) close to Ni 3 5+ 
(to ensure charge neutrality). Using a phenomenological 
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FIG. 2: (color online) (a) Room temperature (300 K) neutron diffraction pattern obtained using HRPD (resolution AQ/Q ~ 
10~~ 3 , 15 hours counting on an 11 g powder sample). The solid curve through data points is a fit to the distorted P6322 space 
group. Vertical bars indicate Bragg peak positions and the bottom curve shows the difference between the fit and the data. 
Zoomed- in regions showing weak supercell peaks are plotted in Fig. [3] (b) Detail of the (111) peak lineshape expected to split 
in the case of a structural distortion to an orthorhombic or monoclinic structure. No splitting could be detected, and the 
lineshape at both 300 K (open symbols) and 2 K (filled circles) could be well described by a resolution-convolved profile for 
the hexagonal P6322 space group (solid lines) with lattice parameters adjusted for thermal contraction upon cooling. 



bond-valence models to relate the valence of the cen- 
tral ion to the bond-lengths, v — J2i e^ r °^ ri ^ B where 
5=0.37 A and r = 1.686 A for the Ni 3 +-0 2 ~ pair, gives 
nominal valences in the ionic limit for the three sites as 
2.42 (Nil) and 3.07 (Ni2 and Ni3), suggesting a signifi- 
cant, but most likely only partial charge disproportiona- 
tion. We note that the bond lengths found in 2H-AgNi02 
are similar to those found in YNi0 3 [d(Nil-0)=1.923 A 
and d(Ni2-0))= 1.994 A, ref. proposed to have a 

charge disproportionation into two Ni sublattices of va- 
lences Ni 3±<T (cr ~ 0.35). 

Upon cooling to lower temperatures no evidence for 
a further structural distortion could be found. At 2 K 
the main and supercell peaks are displaced in wavevec- 
tor following the lattice contraction and the structural 
peaks could be well described [see solid lines in Fig. [3] 
c,f,i)] by the same crystal structure as at 300 K but 
with shorter lattice parameters a — 5.08110(2) A and 
c = 12.24670(7) A (i? B rag g = 4.45%, R F = 5.53%). We 
note that the related triangular-lattice material NaNiC>2 
behaves very differently showing a strong ferro-distortive 
transition into a low-temperature monoclinic structure 
with a significant difference (4%) in the two in-plane lat- 
tice parameters^ We tested for such a scenario in the 
2H-AgNi02 where a departure from hexagonal symmetry 
could be accommodated within the orthorhombic Cmcm 
space group with lattice parameters a x b x c. Within 
the experimental accuracy no splitting of the main peaks 
could be detected [see Fig.O})] and the fitted in-plane lat- 



tice parameters had the same ratio as in the undistorted 
structure (b/a = V3) to within better than 0.02%, so 
we concluded that the hexagonal symmetry is preserved 
down to the lowest temperature probed of 2 K. 



IV. TRANSITION TO THE HIGH-SYMMETRY 
STRUCTURE AT HIGH TEMPERATURES 

We also measured the diffraction pattern at high tem- 
peratures motivated by recent resistivity measurements 12 
reporting a weak anomaly near 365 K, proposed to origi- 
nate from a structural transition. We observed that upon 
heating the triple-cell peaks decreased in intensity and 
could not be observed above T5 = 365(3) K, which co- 
incides with the location of the transport anomaly. Fig. 
[5Ja) shows that triple-cell peaks are absent in the 420 
K data and the main structural peaks are only slightly 
displaced in Q due to lattice expansion. The collected 
diffraction patterns at various temperatures were refined 
in the distorted space group P6322 and the obtained 
temperature-dependence of the oxygen displacement e 
away from the high-symmetry position is plotted in Fig. 
[SJb) inset and shows a rather rapid decrease near Tg. 
The lattice parameters increase smoothly with increasing 
temperature with no clear anomalies near the transition 
apart from possibly a small kink in the c lattice constant 
as linear fits to the data below and above 365 K (solid 
and dashed lines in Fig.[5Jb)) give a slightly smaller slope 
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FIG. 3: Zoomed-in regions of the 300 K neutron diffraction pattern a),b),e),h) showing a number of the weak supercell reflections 
disallowed in the ideal P63 /mmc structure and associated with a tripling of the unit cell in the ab plane. Solid lines are the 
calculated profile for the distorted P6322 structure in Fig.|4]and dashed lines show the estimated local background including the 
resolution tails of nearby main structural peaks (data in some subpanels is shifted vertically by indicated amounts for clarity). 
Paired panels b-c),e-f), h-i) show how the supercell peaks are displaced in Q following the lattice contraction upon cooling (the 
base temperature data has a higher overall background as it was collected in a cryostat). Supercell peaks are not observed in 
the X-ray data in bottom panels d,g,j), consistent with the structural modulation involving mainly displacements of the light 
oxygen ions. 



in the low-temperature phase, consistent with the struc- 
tural distortion in this phase making the crystal lattice 
more rigid so less able to expand upon increasing temper- 
ature. The disappearance of the supercell peaks at high 
temperature indicates a transition to the high-symmetry 
crystal structure where all Ni-0 bonds become equiva- 
lent and charge is uniformly distributed on all Ni sites. 
Indeed the 420 K data could be well described by the 
undistorted Y^jmmc space group with do = 2.94267(3) 
A, c = 12.2554(2) A and oxygen height z Q = 0.07991(6) 
(iWgg = 11.8 %, Rf = 7.13 %). 

To conclude the analysis of the diffraction pattern we 
note that in addition to the triple cell peaks identified 
above, the data also showed diffraction peaks due to a 
small admixture (1%) of the rhombohedral 3R poly type 
(currently at the limit at which pure hexagonal 2H poly- 
type can be chemically prepared) and some other small 
peaks that could not be indexed by any obvious commen- 
surate fractional index of the main peaks and which were 
still present at high temperatures above the structural 
transition at 365 K. Those were attributed to unidentified 
impurity phases below the 1% level that occurred during 
chemical synthesis (scanning electron microscopy (SEM) 
coupled with energy-dispersive X-ray analysis (EDX) of a 



TABLE II: Basis vectors for the irreducible representations 
for charge (scalar) order at the Ni sites (2a) in space group 
PG^/mmc (no. 194) for propagation vector qo — (1/3, 1/3,0) 
obtained using group theory^. 



position 


(0, 0, 0) 


(0, 0, ±) 




1 


1 


r 2 


1 


-1 



small part of the sample indicated small traces of Au, Al 
and Si, but the chemical composition could not be pre- 
cisely identified). The alternative origin of those small 
peaks could have been instrumental (e.g. spurious sam- 
ple holder reflections). 



V. SYMMETRY ANALYSIS OF THE 
STRUCTURAL DISTORTION 

To understand better the mechanism of the struc- 
tural transition at T5 = 365 K from the ideal to 
the distorted structure we performed an analysis of 
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TABLE III: Basis vectors of the irreducible representations for displacement (polar vector) order at the oxygen sites (4/) in 
space group YQs/mmc (no. 194) for propagation vector go = (1/3, 1/3,0) obtained using group theory^. 



position 


(3, 3, zo) 


(|. f > i-zo) 


I- 2 1 1 \ 
1 3 ' 3 ' 2 z °> 
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9 Xfn 
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r^Jj Ni3 
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layer 2 
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V Nil -0 2.022 A 
Ni2-0 1.93* A 
4 Ni3-0 1.934 A 



FIG. 4: (color online) (top) Schematic diagram of the NiCh 
layer at z = \ showing how the displacements (small arrows) 
of the oxygen ions (small balls) lead to a periodic arrangement 
of expanded (large circle, Nil) and contracted (small circles, 
Ni2,3) NiOe octahedra. Thick hexagonal contour shows the 
honeycomb network of contracted sites. The origin of the 
coordinate system is at the circled Ni2 site, (bottom) The ex- 
panded site Nil has a staggered zig-zag arrangement between 
even and odd layers stacked along the c-axis. Layer 2 in the 
unit cell (z = | and — j) is obtained from layer 1 by 180° 
rotation around the central (5,5,2) axis followed by a c/2 
translation. 



symmetry-allowed order patterns in the undistorted, 
high-temperature space group P63 /mmc with the ob- 
served ordering wavevector q — (1/3,1/3,0). We con- 
sidered both the charge (scalar) order on the Ni sites as 
well as the oxygen ion displacements (polar vector order). 
Our refined model for the structure in the distorted 



phase in Fig. 2]shows a periodic arrangement of expanded 
(Nil) Ni 3 ~ CT and contracted (Ni2 and Ni3) Ni 3+0 ' 5,T sites 
with a = 1 in case of complete charge disproportiona- 
tion. To try to reproduce this we looked for symmetry- 
allowed charge (scalar) order patterns with propagation 
vector qo = (1/3,1/3,0) at the Ni 2a sites (2 atoms) 
in the unit cell of the YQ^/mmc space group. The ob- 
tained irreducible representations and basis vectors are 
listed in Table HT1 There are two basis vectors (r^a corre- 
sponding to irreducible representations in Table HI]) 
which physically correspond to having the same or oppo- 
site charges for Ni sites above each other in the two layers. 
The order pattern implied by the structural refinement 
cannot be described by a single basis vector, but can be 
described by a linear combination of two basis vectors, 
i.e. 



c 



J-k/3 



(r 1 +r 2 ) + e l47r / 3 (n - r 2 ) 



(1) 



where Tj(l) is the charge on atom 1 (layer 1) for the basis 
vector Ti etc. 

The corresponding charge order pattern (relative to 
the uniform high-temperature phase where each Ni site 
has valence +3) is shown in Fig. [H] and consists of a tri- 
angular lattice of charges —a, i.e. electron-rich Ni 3_,T , 
surrounded by a honeycombe network of charges +0.5er, 
i.e. Ni 3+0 ' 5<T . This charge order pattern is obtained from 
the basis vector ( modulated by the phase factor due to 
the propagation vector qo, i.e. in the unit cell at distance 
r = niao+n2bo+ n 3 c from the origin (7^2,3 integers) the 
charges are given by the real part of the complex vector 



Qr 



^qar 



For atom 1 (layer 1) the charge is 

Q r (i) = cr cos (qor + 7r/3) 
and for atom 2 (layer 2) is 

Q r (2) = a cos (qor + 5n/3) . 



(2) 



(3) 



(4) 



The normalization prefactor in eq. ((2|) gives the magni- 
tude of the charge order and was chosen such as to obtain 
charge — a on the Nil sites. 

Subsequently we looked at possible order patterns for 
oxygen ion displacements corresponding to the same 
propagation wavevector q . There are in total 4 one- 
dimensional and 2 two-dimensional irreducible represen- 
tations that could describe displacement (polar vector) 
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FIG. 5: (a) Comparison between 300 K (open symbols, lower trace) and 420 K data (filled symbols, upper trace) showing the 
absence of the triple-cell peaks (021), (212) and (314) at high temperature. Solid lines are fits to the distorted (300 K) and 
ideal structure (420 K), respectively and dashed lines show the estimated background level. The high-temperature data has 
been shifted vertically by the indicated offsets for clarity, (b) Temperature dependence of the lattice parameters: c left axis, 
a right axis. Solid (dashed) lines for the c parameter are straight line fits to the data below (above) Ts = 365 K. Filled/open 
symbols are data from different instruments. Inset: temperature-dependence of the oxygen displacement parameter e, solid line 
is guide to the eye. 




layer 2 




FIG. 6: (color online) Charge order pattern described by eq. 
p| to be compared with Fig. [4] Labels -1 and +0.5 indicate 
charges in units of a. Thick hexagonal contours indicate the 
honeycomb network of the electron depleted Ni2,3 sites and 
dashed line contour is the unit cell of the charge-ordered struc- 
ture. Light shaded area is the unit cell of the ideal structure 
with all Ni sites equivalent. 



order at the 4/ oxygen sites (4 atoms per unit cell) and 
the two one-dimensional representations relevant for our 
discussion are listed in Table IIIII Basis vectors of these 
representations transform in the same way under the 
symmetry operations of the PQ^/mmc space group as 
the basis vectors for charge order listed in Table |TT] so we 
use the same symbols t± and r 2 . Here they have complex 
components indicating a degeneracy with respect to rota- 
tion of the displacement vector in the ab plane. We again 
find that the observed structural displacements cannot be 



described by a single basis vector, but is described by a 
linear combination of two basis vectors, i.e. 



1 



,i77T/6 



(ti+t 2 ) 



3 i47t/3 



(ti-t 2 ; 



(5) 



where the complex phase factor in front serves to rotate 
the oxygen displacement vector in the ab plane and with 
the chosen phase the displacement pattern is as in Fig. 
[2J i.e. oxygens displaced radially out of the expanded Ni 
sites (the numerical prefactor is used for normalization). 
In the above equation Ti(l) is the complex displacement 
vector for atom 1 in the representation Ti as per Table 
HTT1 (Note that we use bold symbols to indicate vectors 
with components along the three crystallographic direc- 
tions). The relative phase factor of the t\ and r 2 basis 
vectors is the same for the polar order of oxygen displace- 
ments eq. (O and for the charge order on the Ni sites eq. 
(fTj) as expected if the structural transition was driven by 
charge fluctuations coupled to an oxygen phonon breath- 
ing mode. It is noteworthy that the obtained oxygen 
displacements are consistent with an optic phonon mode 
at the Brillouin zone corner point qa predicted by lattice 
dynamics calculations^, for a structure identical to one 
undistorted Ni0 2 layer. So our results provide support 
that it is this phonon mode (probably modified slightly 
due to couplings with the Ag layer above and below) that 
mediates the V3x \/3 charge order in the triangular Ni0 2 
planes. 

In appendix B we provide further symmetry analysis 
of the basis vector modes for the charge and displace- 
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ment orders using co-representation analysis and prove 
that they uniquely identify the distorted space group as 
P6 3 22. 

To conclude this section for completeness we quote be- 
low the expressions for obtaining the individual atomic 
displacements starting from the basis vector mode The 
order pattern is given by £ modulated by the phase factor 
due to the propagation vector qo, i.e. 

d r = 2e e lqar (6) 

where r = (7^1,^2,%) defines the unit cell and £(1) is 
the basis vector for atom 1 etc. The scale prefactor 2e 
was introduced such as to give the absolute magnitude 
of the final displacements e\/3ao as per Table |TJ The real 
part of the above equation gives the actual displacement 
vectors 

2e (cos (q r + <p x ) , cos (q r + <t>y) , 0) (7) 

where the phases 4> x ^ y are given by the basis vector for 
that atom £ = (e^*,'e**»,0) using eq. © and Table HTT1 
For example for the oxygen atom 1 in the unit cell at the 
origin r = (0, 0, 0) the displacement is 

2e (cos (4tt/3) , cos (5tt/3) , 0) = e(-l, 1, 0). (8) 

VI. SUSCEPTIBILITY AND SPECIFIC HEAT 

Magnetic susceptibility measured on a small powder 
sample is shown in Fig. [3(a) . The pronounced drop at 
low temperatures near 20 K [see inset] is attributed to the 
onset of magnetic order and at higher temperatures above 
~ 70 K the data can be well described by a local-moment 
Curie- Weiss form plus a small temperature-independent 
part xo, 

C 

X = 7f—Q + Xo (9) 

where C = 0.445(5) emu K/mole, 9 CW = -107(2) K, 
and Xo — 1-7(1) X 10~ 4 emu/mole, similar to previous 
reports. 15 The large negative Curie-Weiss temperature 
shows dominant antiferromagnetic interactions. Mag- 
netic order occurs only at significantly lower tempera- 
tures, ~20 K, suggesting that fluctuations due the low- 
dimensionality (mainly in-plane interactions) and the 
frustrated triangular geometry are important in sup- 
pressing the magnetic ordering temperature. 

The effective magnetic moment extracted from the 
Curie- Weiss fit is /x e ff = 1-88 /is per Ni ion, and this was 
used previously^ as evidence that Ni sites were in the 
low-spin Ni 3 + state {tl g e\) with S = 1/2 (^ cff = 1.73// B 
for <? = 2). However, a charge disproportionation sce- 
nario as suggested by our structural measurements could 
also lead to similar values of the effective moment. In the 
extreme case of complete charge disproportionation the 
Nil site is Ni 2+ ^g e g with a large spin moment S = 1, 
whereas the other sites are Ni 3 5+ likely to have only a 



very small spin moment as they are close to Ni 4+ with 
a filled t\ g level, which has 5 = 0. The average effec- 
tive moment observed by high-temperature susceptibility 
considering only S = 1 moments on the Nil sites would 
then be jU e a = 1.63/j.b (for g = 2), which is only 6% lower 
than the value for S = 1/2 at every site, and it is possible 
that including a more realistic scenario of partial charge 
disproportionation (more likely to be the case here) could 
bring this estimate closer to the experiment. 

Finally, we plot in Fig. 0(b) specific heat measurements 
which observe a sharp lambda-like peak identified with 
the magnetic transition near 20 K. No other anomalies 
could be observed up to the highest temperature stud- 
ied of 240 K, consistent with no additional (structural) 
transitions occurring in this temperature range. 



VII. MAGNETIC STRUCTURE 

The magnetic susceptibility shown in Fig. [7] (a) inset 
has a sharp downturn below 20 K as characteristic of a 
transition to antiferromagnetic order and below this tem- 
perature additional reflections are observed in the neu- 
tron diffraction pattern at low Q. The magnetic order 
peaks are most easily seen in the difference pattern 4 K 
— 300 K shown in Fig. [8] and can be indexed with respect 
to the P6322 supercell by the commensurate propagation 
vector k = (1/2,0,0). The magnetic order parameter 
curves are plotted in Fig. [5] as a function of reduced tem- 
perature T/Tn. The Neel temperature obtained from the 
neutron data is 23.7(3) K is slightly higher than the lo- 
cation of the specific heat anomaly near 19.7(3) K, and 
susceptibility maximum near 21(1) K. This is probably 
due to a small temperature calibration offset in the neu- 
tron measurement where the sensor was at some distance 
away from the sample. We regard the absolute value 
of the transition temperature observed by the specific 
heat measurement Tjv=19.7(3) K as the most accurate, 
in agreement with recent /iSR results^ 4 - 

To find the magnetic structure we first consider in Ta- 
ble IIVI the magnetic basis vectors compatible with the 
symmetry of the P6322 crystal structure for the six Ni 
ions in the unit cell (three sublattices) . Symmetry con- 
strains the moments on the same sublattice to be either 
parallel or antiparallcl between the two layers in the unit 
cell. The best fit to the observed diffraction pattern 
(varying the spin direction and magnitude independently 
on the three sublattices) was obtained for the case when 
only one of the three sublattices was ordered, either Nil 
or Ni3, with a moment of 1.552(7) hb along the c-axis and 
spins parallel between adjacent layers (irreducible repre- 
sentation A3 in Table HVl with v 7^ or v" ^ and all the 
other spin components equal to 0). This correlates well 
with the structural analysis whereby the Nil site should 
have a large spin moment (S — 1 at complete charge dis- 
proportionation Ni 2+ ) whereas the other two sites (Ni2 
and Ni3) with valences close to Ni 3 5+ would have a much 
smaller moment, as their configuration would be close to 
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FIG. 7: (a) Inverse magnetic susceptibility (1/x) fitted to a 
Curie- Weiss law eq. (1) (solid line) gives a large negative inter- 
cept indicating dominant antiferromagnetic interactions. In- 
set shows suppression of susceptibility below 20 K attributed 
to onset of antiferromagnetic order, (b) Temperature depen- 
dence of the specific heat observing a sharp lambda-like peak 
near the magnetic transition temperature. 




FIG. 8: Difference pattern 4 K - 300 K from OSIRIS showing 
peaks of magnetic origin, indexed by the propagation vector 
fc=(l/2,0,0). The circles represent the observed intensities, 
the solid curve is a fit to the magnetic structure depicted in 
Fig. [El an d vertical bars indicate the magnetic Bragg peak 
positions. The bottom curve shows the difference between 
the fit and the data. 
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FIG. 9: Observed intensity of several magnetic Bragg peaks 
as a function of reduced temperature T/Tn- Data points are 
from elastic neutron powder scattering measurements using 
IN6. Lines are guides to the eye. 



Ni 4+ with S = 0. Therefore we associate the magnetic 
order with site Nil carrying the largest spin moment. 

The experimental data doesn't preclude the possibility 
of having also a very small magnetic moment along the 
c-axis on either Ni2 or Ni3 sites, these however could 
only be of the order of 0.1 /is at maximum (the spin 
arrangement would then be described by A3 with v' or 
v" equal to —0.1 /is, respectively). Because taking into 
account these small magnetic moments doesn't improve 
the quality of fit in the following we assume for simplicity 
that only the sublattice Nil is ordered with a large spin 
moment. 

The magnetic structure is illustrated in Fig. [TOT a) and 
consists of alternating ferromagnetic stripes, i.e. spins 
are parallel along one of the three directions in the tri- 
angular plane and antiparallel along the other two direc- 
tions. This structure can occur in three different domains 
obtained by ±60° rotation around the c-axis, which 
would correspond to the symmetry-equivalent wavevec- 
tors fc'=(0,l/2,0) and fc"=(l/2,-l/2,0) of the star of k. 
We note that the magnetic structure has only a 2-fold 
symmetry rotation axis along the c-axis whereas the crys- 
tal structure has 3-fold rotation symmetry along the o 
axis. Three equal- weighted domains of those three struc- 
tures would be expected in a macroscopic sample and 
each domain has the same powder-averaged diffraction 
pattern. 

In the following we discuss possible mechanisms to 
stabilize the observed magnetic structure in a model of 
stacked triangular layers, so we consider the Hamiltonian 

n = ^JSi-Sj+J2 J'SrSk+YsS'^-Si-DY,^) 2 ■ 

NN NNN interlayer i 

( 10 ) 

Here NN indicates summing over all in-plane nearest- 
neighbor pairs with coupling J, NNN denotes next- 
nearest-neighbour in-plane pairs with coupling J' [see 
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Fig. flOTb')] and J" is the inter-layer coupling (3 neigh- 
bors above and 3 below, see Fig. [TUta)). with J, J' > 
(antiferromagnetic) and J" < (ferromagnetic). The 
last term in eq. (|10[) is an on site easy-axis anisotropy 
proposed to arise from crystal-field effects and required 
to stabilize the ordering spin direction along the crystal- 
lographic c-axis. 

Considering first the 2D antiferromagnetic Heisenberg 
model (D = J" = in eq. (JTUJ)) the ground state for 
NN couplings only is the 3-sublattice 120° spiral, how- 
ever adding moderate antiferromagnetic NNN couplings 
1/8 < J' I J < 1 stabilizes the collinear stripe order j2£ In 
fact for this model at the classical level the two-sublatticc 
stripe order is degenerate with a continuous manifold of 
4-sublattice non-collinear states but zero-point quantum 
fluctuations are predicted to lift this degeneracy through 
"order by disorder" and select the stripe ground state. 
An easy axis anisotropy D > is expected to further 
stabilize the collinear state thereby reducing the mini- 
mal required J'. At the classical level in the limit of 
very large anisotropy approaching the Ising limit the re- 
quired J' becomes infinitesimally small, but is non zero^ 
In particular the NNN couplings are required to lift the 
degeneracy between the stripe state and the collinear 
state with 2 spins up and one down for each triangle, 
the so-called up-up-down state (UUD), with the ground 
state energy per spin e s = (-J — J')S 2 lower than 
e uuD _ j _j_ 3^)52 for an tiferromagnetic J'. As an 
alternative to NNN couplings we note that the interlayer 
interactions could also provide a mechanism to lift this 
degeneracy as the stacking of magnetic layers energet- 
ically favors the stripe order. The interlayer energy is 
J" S 2 (energy gain for J" < ferromagnetic) for the 
stacked stripe order depicted in Fig. [TUT a) where each 
spin has 4 favorable and 2 unfavorable interlayer bonds, 
whereas for a stacked UUD structure the interlayer en- 
ergy is reduced to J"S 2 /3 because only 2/3 rds of sites 
have 2 net favorable bonds and the remaining l/3 rd have 
2 net unfavorable interlayer bonds. From this we con- 
clude that the observed magnetic structure could be ex- 
plained starting from an easy-axis nearest-neighbor tri- 
angular antiferromagnet with additional weak in-plane 
second neighbor antiferromagnetic exchange or weak fer- 
romagnetic interlayer couplings. 



VIII. DISCUSSION AND CONCLUSIONS 

To summarize, we have reported high-resolution neu- 
tron powder diffraction measurements in the orbitally- 
degenerate frustrated triangular magnet 2H-AgNi02. We 
have observed a set of weak structural reflections, unde- 
tected in previous x-ray measurements, which indicate a 
small structural distortion with a tripling of the unit cell 
in the hexagonal basal plane. We have proposed that 
this could be explained by a periodic contraction and ex- 
pansion of NiOg octahedra in a three-sublattice structure 
as a consequence of charge disproportionation on the Ni 




FIG. 10: (a) Magnetic structure of the Nil sublattice con- 
sisting of alternating rows (dashed lines) of ferromagnetically 
aligned spins. ± symbols indicate the projection of the spin 
moment along the c-axis, dots represent the unordered (Ni2 
and Ni3) Ni sites in the unit cell (solid contour). Thick black 
symbols indicate the moments in the bottom MO2 plane 
[z = j) whereas faint grey symbols correspond to the mo- 
ments in the upper plane (z = 3/4) obtained by shifting the 
pattern by an in-plane offset do- The arrowed line labelled 
J" indicates one of the three inter-layer exchange paths. The 
drawn structure has propagation vector k = (1/2, 0, 0) (Bragg 
peaks indicated by black stars in Fig. 1 1 1 p , equivalent struc- 
tures are obtained by ±60° rotation around the (§,5, z)-axis. 
(b) In a single layer the ordered sites form a triangular lat- 
tice of spacing a. Short and long arrowed lines indicate the 
nearest- and next-nearest-neighbour exchanges J and J' in a 
minimal model proposed to explain the stability of the ob- 
served structure. 




FIG. 11: Reciprocal basal plane showing locations of struc- 
tural (filled circles) and magnetic reflections (stars). Large 
filled circles and large bold hexagon indicate the zone cen- 
ters and Brillouin zone of the ideal, undistorted struc- 
ture (P63 /mmc), whereas small filled circles and thin line 
hexagons are respectively structural supercell peaks and the 
Brillouin zone of the distorted crystal structure (P6322). 
Filled stars are the magnetic Bragg peak positions from 
the magnetic structure with propagation wavevector k — 
(1/2,0,0) shown in Fig. 1101 open stars are peaks from do- 
mains rotated by ±60°. 
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TABLE IV: Basis vectors for irreducible representations for magnetic order at the six Ni ions in the unit cell (3 independent 
sublattices) for a structure with propagation vector fc=(l/2,0,0) for the P6322 space group obtained using group theory (MODY 
package 22 ), u, v, u , v' , u" and v" are independent spin components. 



site 


Nil (2c) 


Ni2 (26) 


Ni3 (2d) 


position 


ri 2 in 

\3> 3 ' 4' 


(- - -) 

V 3 ' 3 ' 4> 


(0, 0, 1) 


(0, 0, |) 


ri 2 In 

V 3 ' 3 ' A> 


(- 1 i) 

V 3 ' 3 ' 4 J 


Ai 


(2m,u,0) 


(2w,m,0) 


(2w',u',0) 


(-2m'-«',0) 


(2u",u",0) 


(2m",«",0) 


A 2 


(0,-u,v) 


(0,-u,-v) 


(Q,-u',v') 


(0,m>') 


(0,-w",«") 


(0,-m" ,-«") 


A 3 


(0,-u,v) 


(0,u,v) 


(Q,-u',v') 


(0,-ti',V) 


(0,-uV) 


(0,«",«") 


A 4 


(2u,w,0) 


(-2u,-u,0) 


(2u',u',0) 


(2u',w',0) 


(2w",u",0) 


(-2u",-w",0) 



sites. We have also observed that the triple cell peaks 
disappear at high temperatures above Ts = 365 K in- 
dicating a structural transition to the ideal undistorted 
structure, where all Ni sites are identical implying that 
charge is uniformly distributed on the Ni sites. The 
low-temperature magnetic diffraction pattern is well ex- 
plained by a structure of ferromagnetic rows ordered an- 
tiferromagnetically, but interestingly with ordered mo- 
ments present only on the electron-rich Ni sites, which 
as a consequence of charge disproportionation would have 
the largest spin moment (5 = 1 if Ni 2+ ). The observed 
magnetic structure on the ordered sites has 2-fold sym- 
metry compared to the 3-fold symmetry of the crystal 
structure. We have proposed that the magnetic struc- 
ture could be explained starting from an easy-axis tri- 
angular lattice antiferromagnet with additional weak in- 
plane next-nearest neighbor antiferromagnetic couplings 
and/or weak ferromagnetic interlayer interactions. De- 
termination of the relative magnitude of the exchanges 
and anisotropy terms requires measurements of the spin 
gap and spin- wave dispersion band- width and such mea- 
surements are in progress^ 

Both the magnetic order and structural distortion ob- 
served here arc very different from the prevailing the- 
oretical model for Jahn- Teller active transition metal 
ions coupled in a triangular lattice arrangement by near 
90° mctal-oxygen-metal bonds, which predicts ferro- 
distortive orbital order and dominant ferromagnetic in- 
plane interactions^ 3 - as indeed observed experimentally 
in NaNi0 2 (Ref. M)- 2H-AgNi0 2 shows a different 
type of structural modulation which suggests an alterna- 
tive mechanism of lifting the large degeneracy in the or- 
bital sector by means of charge ordering, leading to non- 
equivalent Ni sites, some electron rich and others electron 
depleted. In the ideal, undistorted crystal structure each 
Ni site has one electron in the two-fold degenerate e g or- 
bital and at the structural transition below T5 = 365 K 
the lattice separates into two sub-systems. From each 
hexagon of Ni ions one electron jumps to the Ni site in 
the center to form an orbitally-nondegenerate e 2 state, 
and this leaves a surrounding honeycombe network of 
mainly electronically-inactive Ni 4+ sites but with an ex- 
tra electron for every two sites (most likely itinerant and 
therefore distributed with equal probability on every site, 
1/8 filled e g orbital on the honeycombe network). Fur- 



ther experiments in particular NMR, X-ray absorption 
(XAFS) or photoelectron spectroscopy (XPS) would be 
needed to confirm such a scenario and determine quanti- 
tatively the extent of the charge disproportionation. 



Further studies are also needed to understand better 
the very different magnetic behaviour of the expanded 
and contracted Ni sites, attributed in the charge-order 
scenario to electron-rich and -depleted sites, respectively. 
Band structure calculations 16 suggest that as a result of 
charge order the Ni ions in the center of expanded NiOg 
octahedra become more localized and then the magnetic 
order at those sites at low temperatures could be under- 
stood in terms of stacked triangular lattices with a large 
spin moment (valence close to Ni 2+ with 5 = 1), whereas 
the remaining electron-depleted Ni2 and Ni3 sites located 
inside contracted NiOg octahedra maintain a large itin- 
erant character (due to shorter Ni-0 distances and thus 
stronger orbital overlap with the oxygens); in this sce- 
nario the itinerant sites do not show a strong tendency 
to magnetic order because of insufficiently large density 
of states at the Fermi level. In fact, the band structure 
calculations suggested that a small ordered moment of 
~ 0.1 fiB may be induced on the contracted Ni3 sites by 
the ordering of the large Nil moments. We note that 
recent fiSK measurements 2 ^ have reported an anoma- 
lous temperature-dependence of the local static magnetic 
fields which may be due to small ordered moments on 
the honeycombe Ni sites with a different temperature- 
dependence compared to the large moments on the Nil 
sublattice. 



Another interesting aspect worthwhile looking into is 
the role of interlayer coupling in stabilizing the magnetic 
order and structural distortion. It could be investigated 
by looking at a system in which that coupling could be 
modified without modifying the intralayer exchange in- 
tegrals. AgNiC>2 is a rare example of a delafossite having 
two polymorphs whose synthezis is feasible 30 and there- 
fore the ideal candidate for further studies seems to be 
the 3R-AgNi02 polytype. The only difference between 
the latter and the 2H-AgNiC>2 is the way the MO2 layers 
are stacked on top of each other. 
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APPENDIX A: STRUCTURE FACTORS 



- Each of the two irreps generates exactly one corep, 

which has the same matrices as the original irrep on 
the unitary operators. For this reason, we will still 
use the notation t\ and t 2 to indicate the generated 
coreps. 

- Operators hi, h%, h$, /17, hg, hu, Khi^, Khi$, Khyj, 

Khig, Kh 2 i, Kfi23 : all coefficients are 1 for both 
Ti and t 2 

- Operators h 16 , his, hu, h 2 2, h 2 4, h 2n , Kh 4 , Kh 6 , Kh 2 , 

Khio, Khi 2 , Khg: coefficients are 1 for t\ and -1 
for r 2 

The corep modes 8 can be obtained from the irrep 
modes r as 



8 = ^(t + KIt) 



(Bl) 



Here we list the measured and fitted nuclear (Table fV|) 
and magnetic (Table IVTf structure factors. 



Scalar Modes 



APPENDIX B: CO-REPRESENTATION 
ANALYSIS OF THE TRIPLE-CELL CRYSTAL 
STRUCTURE 

Here we look in more detail at the symmetry properties 
of the experimentally determined supercell modulation 
and show how the the basis vectors for the charge and 
displacement order patterns can be used to uniquely iden- 
tify the distorted space group. Specifically we will find all 
symmetry operations that leave the order pattern invari- 
ant to determine the space group. We follow closely ref. 
HH where the basic concepts of co-representation anal- 
ysis as applied to the symmetry reduction for q-vector 
modulations at a generic point in the Brillouin zone are 
explained. We first construct from the irreducible rep- 
resentation (irrep) modes a new set of modes, known 
as corepresentation (corep) modes, which are invariant 
upon application of the anti-unitary operator KI, where 
/ is the inversion operator (/113 in Kovalev notation) and 
K is the complex conjugation. When combined with 
their complex conjugate, these modes are centrosymmet- 
ric even if the propagation vector is not equivalent to 
its inverse, as it is the case here. In other words, the 
"corep little group" contains all 24 symmetry operators 
in PGz/mmc. We can then re- write both the scalar and 
the polar vector modulation as a linear combinations of 
corep modes, and directly assess their symmetry. It is 
important to remember that the coefficients of the corep 
modes are complex-conjugated upon application of an 
anti-unitary operator. 

The irrep-corep matrix for space group 194 and prop- 
agation vector q = (1/3,1/3,0) is listed in Kovalev 29 , 
and, for irreducible representations t\ and t 2 can be sum- 
marized as follows: 



Since Ni is on the inversion center, and the mode is 
scalar and real, we have KIt = r, and 8 = r. We can 
then re-write the scalar modulation in Eq. ^ as 



C = e M/3 \(8i + S 2 ) + e u ^ 3 (Si - S 2 ) 



(B2) 



2. Vector Modes 

Here, site permutation, polar vector inversion and com- 
plex conjugation all come into play, since 



KIt(\) = -- r(2)e i47r/3 
KIt(3) = -7(4)e i4 " /3 



(B3) 



where t(1) is the mode on atom 1 etc. The irrep modes 
are 



t(1) - ±t(3) = (e 



r(2) = ±t(4) = (e 



_ ft,™/! 



t/2 
t/6 



0) 
0) 



(B4) 



where the + and — are for irrep 1 and 2, respectively (see 
Ti and t 2 in Table |TTT|. By inserting eq. JEH} and JB4} 
into eq. (|B1[) , and after some manipulation, one obtains 



V3 
2 



8 = l^ e -™ /& T 



(B5) 



We can then re-write the vector modulation eq. © as 



£ = ^U 4 ^ 3 \{8i + 62) + e^/ 3 (8! - S2) 



(B6) 
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3. Symmetry analysis 



so neither mode is invariant. 



It is now straightforward to perform the symmetry 
analysis on the scalar and vector modulation. There are 
four cases, depending on whether the operators are uni- 
tary or anti-unitary and whether their matrices have the 
same or opposite signs for t\ and r 2 or T\ and r 2 . 

Case 1 : Unitary, same sign (e.g., /i 3 ) 



h 3 ( 



oW3 



(6 1 +6 2 )+e"*' 3 (6 1 -6 2 j\ =C (B7) 



h 3 £ = ^e^ 3 



(S 1 +S 2 ) + e ii ^ 3 {S 1 -S 2 ) 

so both modes are invariant. 
Case 2 : Unitary, opposite sign (e.g., h\e) 



t/3 



(5 1 -5 2 ) + e^/ 3 {5 1 +5 2 )\^C (B8) 



= ^-e^' 3 \(S 1 - 8 2 ) + e i ^ 3 (S 1 + S 2 ) ± i 



so neither mode is invariant. 
Case 3 : Antiunitary, same sign (e.g., Khis) 



Kh 13 ( 



t/3 



(5 1 +5 2 )+ e -^/ 3 (5 1 -5 2 ) ^C(B9) 



(J 1 +5 2 )+e- 44 ^ 3 (<5 1 -<5 2 ) 



Case 4 : Antiunitary, opposite sign (e.g., Kh§) 

KKC = e~™' 3 [(5! - 5 2 ) + e- ii ^ 3 (5 1 + 5 2 )] (BIO) 

= e- t57r / 3 (5 1 +5 2 )+e-^/ 3 (5 1 -5 2 ) 
= e t ^ 3 (5 1 + 5 2 ) + e i5w / 3 (5 1 - 5 2 ) = ( 

"(5i-5 2 ) + e- <47r / 3 (5i+5 2 )" 



-i4tt/3 

3 

V3 
3 

V3 
3 



e -i2V3( 5l + 52 ) + e -i47r/3( 5l _ ^ 
s i4^/3 (5i + ^ + e i2V3( 5l _ 52 )j = £ 

so both modes are invariant. 



In summary, both modes are invariant upon appli- 
cation of the operators hi, /13, /15, /17, /19, /in, if/14, 
Khe, Kh 2 , Kh w , Kh\ 2 , Kh$ and the resulting modu- 
lated structure is invariant by the same operators with- 
out the complex conjugation. These operators represent 
all the proper rotations of space group P6 3 /mmc - in 
other words, the six-fold screw axis and the associated 
orthogonal 2-fold axes. Taking into account the loss of 
translational symmetry due to the propagation vector, 
the resulting space group is P6322 with the unit cell of 
dimensions \/3a x \/3a x c. 
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TABLE V: List of supercell structural peaks associated with 
the tripling of the unit cell in the basal plane with the ob- 
served and calculated unit cell structure factors for the model 
of expanded and contracted Ni-O bonds shown in Fig. 2] The 
observed \F\ 2 is the peak intensity corrected for instrumen- 
tal resolution effects, divided by the peak multiplicity and 
normalized per unit cell of the P6322 group. For complete- 
ness a selection of the nominal peaks of the ideal structure 
(P63 /mmc space group) is also given, and the peak indices 
are given in both space groups. 





Q 


(h,k,l) 


(h,k,l) 


observed 


calculated \F\ 




(A" 1 ) 


P6 3 22 


P63 /mmc 


(10- 14 m) 


(10" 14 m) 




1.026 


002 


002 


6.24 


6.36 




2.052 


004 


004 


7.19 


6.76 


eak 


2.469 


110 


100 


1.14 


0.94 


a 


2.523 


111 


101 


5.89 


5.98 


lomin; 


2.675 


112 


102 


5.86 


6.12 


2.911 


113 


103 


2.79 


2.90 




3.080 


006 


006 


4.49 


4.30 




3.212 


114 


104 


5.83 


5.92 




3.562 


115 


105 


6.40 


6.52 




3.947 


116 


106 


11.36 


11.40 




4.105 


008 


008 


5.69 


5.43 




4.278 


300 


110 


18.07 


16.55 




4.359 


117 


107 


5.34 


5.40 




1.425 


010 


2/3 -1/3 


0.15 


0.16 




1.515 


011 


2/3 -1/3 1 


0.24 


0.24 




1.756 


012 


2/3 -1/3 2 


0.05 


0.09 




2.097 


013 


2/3 -1/3 3 


0.02 


0.02 




2.851 


020 


4/3 -2/3 


0.31 


0.35 




2.898 


021 


4/3 -2/3 1 


0.57 


0.53 




2.936 


015 


2/3 -1/3 5 


0.23 


0.23 




3.031 


022 


4/3 -2/3 2 


0.14 


0.19 




3.393 


016 


2/3 -1/3 6 


0.00 


0.16 




3.514 


024 


4/3 -2/3 4 


0.00 


0.15 




3.807 


211 


4/3 1/3 1 


0.35 


0.22 


- peaks 


3.836 


025 


4/3 -2/3 5 


0.53 


0.50 


3.865 


017 


2/3 -1/3 7 


0.46 


0.26 


cd 



3.910 


212 


4/3 1/3 2 


0.60 


0.62 


CO 

a 


4.074 


213 


4/3 1/3 3 


0.37 


0.42 


3 

CO 


4.105 


026 


4/3 -2/3 6 


0.39 


0.35 




4.295 


214 


4/3 1/3 4 


0.73 


0.66 




4.562 


215 


4/3 1/3 5 


0.30 


0.26 




4.586 


027 


4/3 -2/3 7 


0.49 


0.56 




4.998 


028 


4/3 -2/3 8 


0.35 


0.22 




5.166 


311 


5/3 2/3 1 


0.67 


0.67 




5.209 


217 


4/3 1/3 7 


0.37 


0.19 




5.242 


312 


5/3 2/3 2 


0.71 


0.74 




5.366 


313 


5/3 2/3 3 


0.56 


0.48 




5.535 


314 


5/3 2/3 4 


0.77 


0.77 




5.575 


218 


4/3 1/3 8 


0.63 


0.58 




5.703 


040 


8/3 -4/3 


0.54 


0.60 




5.726 


041 


8/3 -4/3 1 


0.89 


0.90 




5.745 


315 


5/3 2/3 5 


0.56 


0.65 




5.794 


042 


8/3 -4/3 2 


0.41 


0.32 




5.821 


1 11 


2/3 -1/3 11 


0.22 


0.21 




5.870 


2 10 


4/3 -2/3 10 


0.17 


0.12 



10 



TABLE VI: List of magnetic Bragg peaks with the observed 
and calculated unit cell structure factors for the magnetic 
structure in Fig. 1101 The observed \F\ 2 is the peak intensity 
corrected for instrumental resolution effects, divided by the 
peak multiplicity and normalized per unit cell of the P6322 
group. For completeness the peak indices are given both in the 
distorted (P6322) and the ideal (PG^/mmc) crystal structures. 



Q 


(a, fe, I) 


(h,k,l) 


„ 1 1 I771I 

observed \r \ 


„„l„..l„j „J I T71| 

calculated \b \ 


( A ) 


r 63ZZ 


P63 / mmc 


/in— 14 m \ 
(1U mj 


(1U mj 


0.714 


1/2 


1/6 1/6 


0.61 


0.71 


0.879 


1/2 1 


1/6 1/6 1 


0.28 


0.33 


1.236 


1/2 -1 


1/2 


0.00 


0.00 


1.250 


1/2 2 


1/6 1/6 2 


0.37 


0.39 


1.339 


1/2 -1 1 


1/2 1 


0.74 


0.72 


1.607 


1/2 -1 2 


1/2 2 


0.00 


0.00 


1.697 


1/2 3 


1/6 1/6 3 


0.16 


0.16 


1.889 


1/2 1 


-1/6 5/6 


0.58 


0.63 


1.957 


1/2 1 1 


1/6 5/6 1 


0.34 


0.35 


1.957 


3/2 -1 1 


1/6 5/6 1 


0.34 


0.35 


1.974 


1/2 -1 3 


1/2 3 


0.48 


0.45 


2.142 


3/2 


1/2 -1 


0.00 


0.00 


2.150 


1/2 1 2 


1/6 5/6 2 


0.50 


0.53 


2.150 


3/2 -1 2 


1/6 5/6 2 


0.50 


0.53 


2.173 


1/2 4 


1/6 1/6 4 


0.15 


0.20 


2.202 


3/2 1 


1/2 -1 1 


0.70 


0.68 



